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ABSTRACT 

In this paper, we present both timing and spectral analysis of the outburst of 4U 
0115+63 in April - May 2008 with INTEGRAL and RXTE observations. We have 
determined the spin period of the neutron star at ^ 3.61430 ± 0.00003 s, and a spin 
up rate during the outburst of P = (—7.24 ± 0.03) x lO^^s d~^, the angle of periapsis 
w = 48.67° ±0.04° in 2008 and its variation (apsidal motion) w = 0.048° ± 0.003°yr-i 
when comparing with previous measured values of ui. We also confirm the relation 
of spin- up torque versus luminosity in this source during the giant outburst. The 
hard X-ray spectral properties of 4U 0115+63 during the outburst are studied with 
INTEGRAL and RXTE. Four cyclotron absorption lines (fundamental at ~ 10 — 15 
keV, 2nd harmonic at ~ 22 keV, 3rd at ^ 34 keV and 4th at ^ 45 keV) are detected 
using the spectra from combined data of IBIS and JEM-X aboard INTEGRAL in the 
energy range of 3 - 100 keV. Two cyclotron absorption lines are measured from the 
spectra combining PGA and HEXTE data in the energy band of 4 - 50 keV. The 5-50 
keV luminosities at an assumed distance of 7 kpc are determined to be in the range of 
(1.5 — 12) X lO^^ergs during the outburst. The fundamental absorption line energy 
varies during the outburst: around 15 keV during the rising phase, and transiting to 
~ 10 keV during the peak of the outburst, and further coming back to ~ 15 keV 
during the decreasing phase. The variations of photon index show the correlation with 
the fundamental line energy changes: the source becomes harder around the peak of 
the outburst and softer in both rising and decreasing phases. This correlation and 
transition processes during the outburst need further studies in both observations 
and theoretical work. The known relation of the fundamental line energy and X-ray 
luminosity is confirmed by our results, however, our discoveries suggest that some 
other factors besides luminosity play the important role in fundamental line energy 
variations and spectral transitions. 

Key words: pulsars: individual (4U 0115+63) ~ stars: neutron ~ stars : binaries : 
close - X-rays: binaries. 



1 INTRODUCTION 

The X-ray transient 4U 0115+63 is one of the best studied 
X-ray sources with different instruments. It is a high mass 
X-ray binary composed of a neutron star and a Be star. 
Ever since it was first observed in the Uhuru satellite sur- 
vey (Giacconni et al. 1972, Forman et al. 1978), series of 
observations have been done to study its different kinds of 
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features. With the precise positional determinations by the 
satellites SAS-3, Ariel V and HEAO-1, the companion star 
was identified as a heavily reddened 09e star, V635 Cas, 
with visual magnitude V « 15.5 (see Cominsky et al. 1978, 
Johnston et al. 1987, Johns et al. 1978, Hutchings & Cramp- 
ton 1981). The Be star is an early type, non-supergiant star 
which some times exhibits Balmer series lines. The distance 
of the star is suggested to be ~ 7 kpc (Negueruela & Okazaki 
2001). The spin period of neutron star in 4U 0115+63 was 
discovered at ~ 3.6 s (Cominsky et al. 1978). The binary's 
orbital elements were first determined by Cominsky et al. 
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(1978), Rappaport et al. (1978), and Kelly et al. (1981), us- 
ing the SAS — 3 data sets, finding the eccentric(e = 0.34), 
and orbital period 24.3 days, and precise timing analysis 
in different time intervals have given a upper limit on the 
apsidal motion Co < 2.1°yr^^. 

This hard X-ray emission is produced via accretion pro- 
cess by a neutron star accreting material from its companion 
Be star. However, due to different physical and geometrical 
conditions in the accretion disc, the source has shown dif- 
ferent types of X-ray activities. In quiescence, 4U 0115-1-63 
exhibits low luminosity (La, < 10^* ergs s~^) (Campana et al. 
2001). Some times the object presents Type I X-ray out- 
bursts (lO^^ergs s~^ < Lx < lO^^ergs s~^) separated by 
orbital period. For this source, the type II outbursts ap- 
pear more often then the Type 1. It shows giant X-ray 
outburst(I/2, > lO^'^ergs s~^) enduring for several weeks or 
months. The onset of the Type II outbursts occur near peri- 
astron, and some times it is understood as quasi-periodic be- 
haviors with a modulation period of about 3 to 5 years. The 
physical origin of the outburst mechanism is still unclear. 
Negueruela & Okazaki (2001) suggest that the distance at 
which the circumstellar disc is truncated depends on the or- 
bital elements and viscosity. Binary system with low eccen- 
tricity would expect the gap between the disc outer radius 
and the lobe radius of the Be star so large that normally the 
neutron star could not accrete enough gas at periastron pas- 
sage to express Type I outburst. In this case, these systems 
would only display occasional Type II outburst. 

Research about cyclotron resonance spectral features 
(CRSFs) of this source has been carried out repeatedly with 
different observations since CRSF was clearly detected as 
an absorption feature at ~ 20 kcV (Whcaton 1979) using 
the data from HEAO-1 A4. Soon after that, White et al. 
(1983) pointed out that the CRSF at ~ 20 keV in fact was 
the second harmonic resonance, and two cyclotron absorp- 
tion lines at the fundamental energy ~ 11.5 kcV and 23 
keV were detected. The double harmonic absorption fea- 
tures were confirmed using Ginga data (Nagase et al. 1991; 
Tamura et al. 1992). Heindl et al. (1999) discovered a third 
absorption line at ~ 33.56 keV in 1999 with RXTE. More- 
over, a fourth absorption line was detected by Santangelo 
et al. (1999) at 49.5 keV with BeppoSAX and Tsygankov 
et al. (2007) at 44.93 keV with INTEGRAL. Recently the 
fifth harmonic resonance absorption line at ~ 53 keV was re- 
ported by RTXE observations (Ferrigno et al. 2009). So that, 
4U 0115-f 63 is the only object in the spectrum of which five 
cyclotron line harmonics were detected, making it as one of 
the best objects for studying the physics of cyclotron reso- 
nance in the polar cap regions of X-ray pulsars in binaries. 

4U 0115-1-63 was observed at different X-ray luminos- 
ity levels, and the fundamental energy of CRSFs was re- 
peatedly measured axound 11 keV around the peak of the 
Type II outbursts. But the variations of the fundamental 
line energy were detected when we observed the source in 
the lower X-ray luminosity. The report from Ginga (Mihara 
1995; Mihara et al. 1998) found that during the outburst 
in 1991, when the X-ray luminosity is about 7 times lower 
than the typical luminosity of the outburst peaks, instead of 
the familiar double absorption features at 11 and 22 kcV, 
the outburst exhibited a single deep and wide absorption at 
~ 16 keV. The further studies suggested that the 11 keV 
fundamental energy line moves to ~ 16 keV in lower lumi- 



nosity ranges instead of that the 22 keV second harmonic 
decreasing to 16 kcV (Nakajima et al. 2006). Detailed studies 
showed a luminosity-dependent change of the fundamental 
absorption line energy (Burnard et al. 1991; Nakajima et al. 
2006; Tsygankov et al. 2007). 

4U 0115-1-63 undergone the giant Type II outbursts fre- 
quently about every 3-5 years. The latest giant outbursts 
occurred in in April - May 2008. During this giant out- 
burst, both INTEGRAL and RXTE have carried out point- 
ing observations on the source frequently. In this paper, we 
analyzed both INTEGRAL and RXTE data of 4U 0115-1-63 
obtained during giant outbursts and studied the temporal 
and spectral properties and variations during the outbursts. 
We will concentrate on the spin properties of the neutron 
star, deriving the orbital elements of the binary and the res- 
onance absorption energy variations during the outburst. 



2 OBSERVATIONS 

The giant outburst of 4U 0115-1-63 occurring in 2008 was ob- 
served by INTEGRAL(JJVTEernationaJ Gamma-Ray Astro- 
physics Laboratory, Winkler et al. 2003) and RXTE(Rossi 
X-ray Timing Explorer, Bradt et al. 1993). 

INTEGRAL is ESA's currently operational space- 
based hard X-ray/soft gamma-ray telescope. We mainly 
use the data collected with the softer gamma-ray imager 
(IBIS/ISGRI; Lcbrun et al. 2003) and the small X-ray tele- 
scope JEM-X (Lund et al. 2003). IBIS has a 12'(FWHM) 
angular resolution and arcmin source location accuracy in 
range of energy at 18 - 200 keV, JEM-X obtains a lower 
energy band at 3 - 35 keV. Two instruments are co-aligned, 
allowing simultaneous observations in a wide energy range. 

RXTE is a space X-ray telescope run by NASA. The 
pointing observations are performed mainly by two instru- 
ments: Proportional Counter Array (PCA) and High Energy 
X-ray Timing Experiment (HEXTE). PCA is a spectrome- 
ter with an effective area of 6400 cm^, with energy range of 
2-60 keV and time resolution of 1 microsec. HEXTE covers 
a higher energy range of 15 - 200 keV with a timing resolu- 
tion of 8 microsec. The two instruments aboard RXTE are 
very effective for timing analysis on fast variability of X-ray 
sources. 

In Table 1, we summarize the archival data of 4U 
0115-F63 during the outburst in 2008 with both INTEGRAL 
and RXTE observations, and the instrument, date and ex- 
posure time are listed in details. 

The used INTEGRAL data are available from the IN- 
TEGRAL Science Data Center (ISDC). The analysis was 
done with the standard INTEGRAL off-line scientific anal- 
ysis (OSA, Goldwurn et al. 2003) software, ver. 9.0. IN- 
TEGRAL observed the outburst in seven revolutions (one 
INTEGRAL orbit, 3 days). Each revolution includes several 
individual pointing observations (one science window) last- 
ing from about 3000 s to 10000 s. Individual pointings in 
each satellite revolution processed with OSA 9.0 were mo- 
saicked to create sky images for the source detection. For 
each revolution, we got one spectrum from 3 - 100 keV us- 
ing the combined spectra of JEM-X and IBIS. IBIS has a 
poor time resolution which is larger than 1 second, not suit- 
able for timing analysis of the 3.6 sec X-ray pulsar in 4U 
0115-1-63. In timing analysis, only JEM-X data are used to 
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Table 1. Observations of 4U 0115+63 during the outburst in 2008 
with INTEGRAL and RXTE, the instrument, observed date, ex- 
posure time are shown. 



Instrument 


MJD (Rev.) 


Exposure time (s) 


JEM-X 


54546.52(664) 


99599.08 


JEM-X 


54555.46(667) 


49290.03 


JEM-X 


54558.45(668) 


39747.02 


JEM-X 


54563.38(669) 


50285.74 


JEM-X 


54566.41(670) 


46669.91 


JEM-X 


54573.90(673) 


52292.43 


JEM-X 


54580.62(675) 


110950.9 


IBIS 


54546.52(664) 


99599.08 


IBIS 


54555.46(667) 


50067.22 


IBIS 


54558.45(668) 


51651.26 


IBIS 


54563.38(669) 


53380.74 


IBIS 


54566.41(670) 


49267.91 


IBIS 


54573.90(673) 


52293.74 


IBIS 


54580.62(675) 


110952.9 


Instrument 


} : : ^ 

MJD (pointing) 


Exposure time (s) 


HEXTE 


54545. 265 (930320 10 100) 


13984 


HEXTE 


54545.470(93032010101) 


1680 


HEXTE 


54548.074(93032010200) 


16928 


HEXTE 


54549.186(93032010201) 


11424 


HEXTE 


54549.971(93032010202) 


14864 


HEXTE 


54550.236(93032010203) 


4256 


HEXTE 


54546.970(93032010204) 


20640 


HEXTE 


54551.083(93032010205) 


7584 


HEXTE 


54551.284(93032010206) 


2336 


TT TP V T^TT^ 


54556.254(93032010300) 


9936 


xlJliA. i H/ 


04000. 14:Z(^yoUozUlUoUl ) 


3088 


nUiA. 1 hj 


o4oo4..^2bt9303.£0i03U.^ ) 


13584 


rlJiA. i ej 


o4ooo.yo.^^youo.^uiuouo ) 


15392 


U "V "C 


o4oo0.04o(yoOo.i0104UO ) 


16096 


XJT7' YT'TT" 

rlUjA. i III 


o4oD.^.UlH_yoUo.^UlU4Ui ) 


14816 


rlrLA. 1 Hi 


o4ooz. y4o (VJoUozU iU4Uz ) 


12944 




54563.986(93032010403) 


8384 


HEXTE 


54564.817(93032010404) 


15504 




o4ooo.oo4^y JUo JU1U4UO ) 


21344 




o4oDD.y ( o^ydUo Juiu4Ub ) 


2208 




o4oD ( .UUHydUd.^UlUoUU ) 


1376 




o4oo7.yoD^yd0d.£0i0o0i ) 


16160 




o4oDy.yod(ydOd.aiDOioo ) 


15520 


TT TP V T^lIP 

hit JV i ±L 


o4o ( i.y44(ydOd.^lD010iJ 


15616 


HEXTE 


54574.037(93032160200) 


9680 


XJT7' Yn^T7 


o4o ( o. ooU (^UoUoz i()UzUi ) 


13136 


PCA 


54545.265(93032010100) 


13984 


PCA 


54545.470(93032010101) 


1680 


PCA 


54548.074(93032010200) 


16928 


PCA 


54549. 186(93032010201) 


11424 


PCA 


54549.971(93032010202) 


14864 


PCA 


54550.236(93032010203) 


4256 


PCA 


54546.970(93032010204) 


20640 


PCA 


54551.083(93032010205) 


7584 


PCA 






PCA 


54556.254(93032010300) 


9936 


PCA 


54556. 142(93032010301) 


3088 


PCA 


54554.226(93032010302) 


13584 


PCA 


54558.932(93032010303) 


15392 


PCA 


54560.043(93032010400) 


16096 


PCA 


54562.011(93032010401) 


14816 


PCA 


54562.945(93032010402) 


12944 


PCA 


54563.986(93032010403) 


8384 


PCA 


54564.817(93032010404) 


15504 


PCA 


54565.864(93032010405) 


21344 


PCA 


54566.975(93032010406) 


2208 


PCA 


54567.001(93032010500) 


1376 


PCA 


54567.956(93032010501) 


16160 


PCA 


54569.983(93032160100) 


16620 


PCA 


54571.944(93032160101) 


15616 


PCA 


54574.037(93032160200) 


9680 


PCA 


54578.360(93032160201) 


13136 



search for spin period of the neutron star in 4U 0115+63. 
The light curves of JEM-X have a time resolution of 0.5 s. 

Data of RXTE arc available from the standard products 
of the RXTE data are available on RXTE archive, which are 
used for the light curve extractions with a time resolution 
of 0.125 s from both of PCA and HEXTE. Each science 



window records data with the durations of several thousand 
seconds (see Table 1). The spectrum for each science window 
is extracted using the combined data of PCA and HEXTE. 



3 DATA ANALYSIS AND RESULTS 

3.1 TIMING ANALYSIS AND ORBITAL 
ELEMENTS OF 4U 0115-1-63 

The spin properties of the accretion pulsars are crucial for 
us to understand the accretion physics and giant outburst 
mechanisms. The determination of spin period and period 
derivatives of the X-ray pulsars in binaries will sensitively 
depend on the orbital parameters. In addition pulse profile 
variations of X-ray pulsars during the outbursts is not well 
understood and requires more analysis. In this subsection, 
we will study the temporal properties of 4U 0115+63 during 
the outburst of 2008 using the INTEGRAL and RXTE data. 
We will try to derive the spin period of the X-ray pulsar 
combined with the orbital parameters. Study the apsidal 
motion of X-ray binaries is one of the important features to 
understand the nature of the X-ray binaries' properties, so 
we also try to search for the apsidal motion of 4U 0115+63 
by comparing the results in different time intervals. 

X-ray Ught curves of 4U 0115+63 during the 2008 out- 
burst by JEM-X are extracted in two energy bands: 3- 10 
keV and 10 - 35 keV, with a time-resolution of 0.5 s. The 
light curves of PCA are used in the energy band of 2 — 40 
keV. 

A general way to perform timing analysis in determin- 
ing the pulse period of X-ray pulsars in binaries is to find 
the pulse arrival time (see Rappaport et al. 1978, Tamura et 
al. 1992). The pulse arrival time for each data set was deter- 
mined by comparing the folded pulse profile with a reference 
profile. This method can get the exact arrival time of each 
pulse of each data set. However, since the pulse profile of 
X-ray pulsar would change during the outburst, the deter- 
mined arrival time of the position of pulse profile (e.g., peak) 
may have large uncertainties, it needs additional statistical 
test to ensure the results. In this paper, we have carries out 
another way to do the timing analysis on pulse period stud- 
ies (see Raichur & Paul 2010). It would not be affected by 
the variation of the pulse profile. 

This method is mainly based on the HEAsoft, published 
by NASA which aims to handle the high energy data set ob- 
tained from the space satellite, such as INTEGRAL, RXTE. 
First of all, we corrected the times of each photon to the 
barycenter of the solar system, then we used EFSEARCH (a 
build-in function in HEAsoft) to calculate the estimated spin 
period of the source in each data set with best value based 
on Chi-Square estimation. The EFSEARCH might depend 
on the epoch we used, so we could correct it as follows: calcu- 
late the with different epochs, T(epoch) — to + 0.1 i P/N, 
where i = 0, 1, 2, 9, and N is the number of phase bins in 
pulse profile. Then we average ten sets of versus period, 
fit the peak with quadratic and find the maximum which 
should be the estimated period. The la error could be given 
when the quadratic is lowered by 1.0. 

As we get all the pulse periods of data sets during the 
outburst, we fit the data according to the doppler motion of 
this source. The observed pulse period is given as below: 
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i(io) + Pspin{t — to)) 



1 + Vr/C 
1 — Vr/c 



(1) 



Here Pgpin denotes the derivative of spin period, P°p/„ is 
the spin period values observed at time t and to is the ref- 
erence time of the outburst. Vr is the component of velocity 
along the sight of neutron star which depends on the orbital 
elements: 



Vr = 



2Tvax sini 



(1 



-(cos(i/ + uj) + ecosoj), 



(2) 



the aa:Sini denotes for the projected semi-major axis, Port 
is the orbital period, uj is the angle of periastron and e is the 
eccentricity of the neutron star's orbit, v is the true anomaly 
and is closely related to the periastron passage time TL,, the 
relationship could be given as follows: 



tan ■ 



tan ■ 



E /l+_e 



esmiJ = — — (t 



(3) 



(4) 



Here E stands for the eccentric anomaly. Combining all 
those equations, we could get the doppler motion of the 
neutron star in the different orbital phases by fitting the 
observed spin period values in different time intervals. 

Using the EFSEARCH tool, we derive the observed val- 
ues of spin period by INTEGRAL/ JEM-X and RXTE/PCA 
from MJD 54545 - 54580 in the 2008 outburst. The ob- 
served data points of the spin period during the outburst 
are presented in Fig. 1. These data points are fitted by 
the function described by Eq. (1). The fitting is based on 
Levenberg-Marquardt algorithm for nonlinear least square 
method, weighted by error of each point. The parame- 
ters' errors are calculated from the covariance matrix of 
the parameter estimates. We adopted some orbital param- 
eters from Bildsten et al. (1997), i.e., the orbital period 
Port. = 24.317037 ± 0.000062 day, a^sini = 140.13 ± 0.08 
It-sec. Then, we finally get the spin period, the derivative 
and other orbital elements of 4U 0115-f63 as shown in Ta- 
ble 2. 

The spin period of the neutron star is determined at 
~ 3.61430 ± 0.00003 s at the day MJD 54566, with the 
derivative of P ~ (-7.24 ± 0.03) x lO"** s day"\ and T„ = 
MJD 54531.7709 ±0.0603. The neutron star in 4U 0115-f63 
still undergone the spin-up process during the 2008 giant 
outburst, which is consistent with the previous results on 
the history outbursts. Here we summarize previous results 
on the spin period and period derivative of 4U 0115-f63 
in Table 3 for a comparison. From 1981 - 2008, the de- 
rived spin period values show the long-term spin-up trend. 
During the history outbursts, during the outburst in 2008, 
P /P ~ —7.31 X 10"** yr~^, which is considerably larger than 
before. 

We also try to measure the apsidal motion of the bi- 
nary system. Using the INTEGRAL and RXTE data in 
the 2008 outburst, we have derived the angle of periastron 
w ~ 48.67° ± 0.04° in 2008. We collected the previous mea- 
surements of the angle of periastron of 4U 0115-1-63 in differ- 
ent epoches in Table 4. In Fig. 2, the angle of periastron uj 
versus the observed date is plotted. The angle of periastron 
seems to change in the last 30 years, suggesting existence of 
the apsidal motion. A linear function is used to fit the data 



Table 2. Spin period and orbital elements of 4U 0115-(-63 at MJD 
54566 



Parameter 


Value 


P (s) (at MJD 54566) 


3.61430 ± 0.00003 


Papin {sday-'^) 


(-7.24 ±0.03) X 10-^ 


^(°) 


48.67 ±0.04 


u](°yr~^) 


0.048 ± 0.003 


{MJD) 


54531.7709 ±0.0603 


Porb(d) 


24.317037 ± 0.000062 


e(l) 


0.3402 ± 0.0002 


ax sini {It — sec)'-') 


140.13 ± 0.08 


Reduced x^{d.o.f.) 


0.85 (68) 


By referring tire results 


in Bildstcn et al. (1997) 


3.617 , , , , , , , , , ,- 





3.613 
0.0002 

15 

= 0.0000 

a. -0.0002 




54540 54545 54550 54555 54560 54565 54570 54575 54580 54585 54590 

MJD(day) 



Figure 1. Spin period variations due to doppler motion 
of neutron star in 4U 0115+63 during the 2008 outburst. 
The black squares denote for the spin period measured by 
INTEGRAL/ JEM-X, and the solid squares by RXTE/PCA. The 
residuals are shown in the lower panel. The solid curve represents 
the expected doppler motion of the binary due to the best fitted 
orbital elements. 



points, also using Levenberg-Marquardt method and derive 
a apsidal motion rate Co = 0.048° ±0.003° yr^^, error is cal- 
culated from the covariance matrix of the linear fit, of 
this fit is 0.86. 

In addition, we also have studied the evolution of pulse 
profile of 4U 0115±63 during this giant outburst. In Figure. 
3, we displayed the pulse profiles with the outburst time in 
two energy bands of JEM-X: 3-10 keV and 10 -35 keV. 
There exist obvious different features in pulse profiles of the 
high energy band and low energy band. In the lower band, 
pulse profiles show secondary peaks in most time of the gi- 
ant outburst, and may disappear near the end of the out- 
burst (MJD 54581.7). As the fiux arrives at its peak, the 
secondary peak's amplitude also increases, so that we would 
see the secondary peak clearly (see the upper panel of Fig- 
ure 3, MJD 54558.6). As the outburst flux decreases, the 
secondary peak also starts to merge with the main peak. 
While at high energy band, the second peak is not so obvi- 
ous. After the outburst reaches its peak, the secondary peak 
could be seen (see lower panel of Figure 3, MJD 54574.0). 
This difference may be due to the different mechanisms of 
pulsation at different energy levels. 

In Fig. 4, we also presented the pulse proflles with the 
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Figure 3. The evolution of pulse profiles of 4U 0115+63 during the outbursts observed by INTEGRAL/ JEM-X in two energy bands: 
3-10 keV and 10-35 keV, respectively 
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Table 3. History of spin period and its derivative of 4U 0115+63 



MJD 




P/P (yr-i) 


References 


40963 


3.614658 


± 0.000036 


(-3.3 ± 1.4) X 10-** 


Kelley et al. (1981) 


42283 


3.6142 


± 0.0001 


-3.4 X 10-5 


Whitlock et al. (1989) 


43540 


3.6145737 


± 0.0000009 


(-3.2 ±0.8) X 10-5 


Rappaport ct al. (1978) 


44589 


3.6146643 


± 0.0000018 


(-2.6 ± 0.5) X 10-" 


Rickctts ct al. (1981) 


47941 


3.614690 


± 0.000002 


(-1.8 ± 0.2) X 10-* 


Tamura et al. (1992) 


49481 


3.6145107 


± 0.0000010 


(-1.8 ±0.6) X 10-* 


Scott et al. (1994) 


50042 


3.614499 


± 0.000004 


(-4.2 ±0.9) X 10-5 


Finger et al. (1995) 


50307 


3.614451 


± 0.000006 




Scott et al. (1996) 


51232 


3.614523 


± 0.000003 




Wilson et al. (1999) 


51240 


3.61447 


± 0.00002 


(-2.52 ±0.6) X 10-4 


Raichur & Paul (2010) 


53254 


3.61436 


± 0.00002 


(-3.84 ±0.55) X 10-* 


Raichur & Paul (2010) 


54566 


3.61430 


± 0.00003 


(-7.31 ±0.03) X 10-" 


this work 



By directly comparing previous spin period in 1971 and 1974 



Table 4. Epoch and oj history of 4U 0115+63 



Orbit number Periastron time passage (TJD) (°) Satellite References 






40963.08 ±0.17 


51.10 ± 3.60 


Uhuru 


Kelly ct al. (1981b) 


106 


43540.451 ±0.006 


47.66 ± 0.17 


SAS-3 


Rappaport ct al. (1978) 


149 


44585.700 (*) 


47.15 ± 0.13 


Ariel-6 


Ricketts et al. (1981) 


287 


47941.530 ±0.006 


48.02 ± 0.11 


Ginga 


Tamura et al. (1992) 


342 


49279.2677± 0.0034 


47.66 ± 0.09 


BATSE 


Cominsky et al. (1994) 


422 


51224.6465± 0.051 


48.50 ± 0.92 


RXTE 


Raichur & Paul (2010) 


505 


53243.038± 0.051 


50.07 ± 1.86 


RXTE 


Raichur & Paul (2010) 


558 


54531. 7709± 0.0603 


48.67 ± 0.04 


INTEGRAL&RXTE 


this work 



^'^^ The orbit number and is derived from Ricketts et al. (1981) 




Figure 4. The evolution of pulse profiles of 4U 0115+63 observed by RXTE/PCA during the outburst in the energy band of 2 - 40 keV. 



Cyclotron Resonance Energies and Orbital Elements of 4U 0115+63 7 



55 
54 
53 
52 

TJ 

49 
48 
47 
46 





1 ' 1 ' 1 




1 ' 




Uhuru 












RXTE 








RXTE 






1 


SAS-3 

5 _ — 

Ariel-6 

I.I.I 


Ginga i 

— 
— """ n 

BATSE 


1 


RXTE& 
INTEGRAL ■ 

1 



40000 42000 44000 46000 480Q0 50000 52000 54000 56000 
Time (MJD) 



Figure 2. The measured u) values of 4U 0115+63 from 1971 to 
2008. A linear function is used to fit the data points in deriving 
average o) ~ 0.048° ± 0.003°. 



outburst evolution by RXTE/PCA in the energy band of 2 
- 40 keV. PCA has a better time resolution (0.125 s) and 
wider energy band, that could show us a vivid picture of the 
evolution of pulse profile during the outburst. The second 
peak was dotoctod in all time intervals, and it gradually 
grows up, and becomes more obvious around the peak of 
the outburst, then it begins to fade, which is consistent with 
the pulse profile evolution detected by INTEGRAL/JEM- 
X. The last pulse profile in MJD 54578.4 (Figure 4) almost 
merges with the main peak. 



3.2 Spectral Analysis and Cyclotron Absorption 
Features 

In this section we performed the spectral analysis of the hard 
X-ray spectra from several kcV up to 100 kcV with RXTE 
and INTEGRAL respectively. Wc ignored spectral data with 
poor signal to noise ratio in several energy channels, and 
used the spectral data as follows: (a) JEM-X: 3-31 keV; (b) 
IBIS: 28-100 kcV; (c) PCA: 4-30 keV; and (d) HEXTE: 17 
- 50 keV. 

We fit the data using the XSPEC package. All the spec- 
tra have been fitted with a power-law type spectra modified 
by a high energy cut off structure (the powerlaw*highecut 
model in XSPEC package; White et al. 1983), additionally 
modified by cyclotron absorption lines in Lorcnz profile (see 

{WjE/E,y^lf ' 
(B-B,„,,)2 + W|_ 

iron line (gaussian profile) at 6.4 keV (see Tsygankov et al. 
2007). Besides, photo-electric absorption (see Morrison and 
McCammon 1983) was added when fitting the INTEGRAL 
spectral data. Hard X-ray flux is measured in the range of 
4-50 keV for both satellites. 

For the spectra combined with JEM-X and IBIS, around 
the peak of the outburst, the lower band (3-4 keV) of spec- 
trum is a little lower than the contirmum model spectrum if 
we only use the power law model with high energy cut off 
to fit the continuum spectra. Hence, we add the column ab- 
sorption component in the fittings. But for the PCA's data. 



Mihara et al. 1990): exp 



-D, 



with a fixed 



because we ignore the energy channels below 4 keV, where 
photo-electric absorption is not prominent in the low energy 
band, therefore column density absorption is not included. 

The spectra of 4U 01154-63 during the outburst ob- 
tained by JEM-X and IBIS are shown in Figure 5. The fit- 
ted spectral parameters are presented in Table 5. For Rev 
675, the detection significance level of the source by JEM- 
X and IBIS is low (< 7cr), no line features could be found 
in the spectrum, so we do not include it in the following 
analysis. For the spectrum in the range of 3 100 kcV by 
INTEGRAL, We observed cyclotron resonance absorption 
features up to the 5th harmonic resonance energy. The fun- 
damental line energy is determined at 10 15 kcV in differ- 
ent cpoches of the outburst. The second harmonic energy is 
around 22 keV, and the third harmonic energy is about 34 
keV, and the fourth harmonic energy is around 42 keV. The 
fifth harmonic energy has a large uncertainty which varies 
from 49 keV to 55.5 keV, but the F-test shows this com- 
ponent might not be so significant (probabilities are rrmch 
higher than 0.01), so we do not include the 5th absorption 
line in the spectral fittings. We also listed the probability of 
F-test, The result obtained through PCA and HEXTE are 
presented in Table 6 & Fig. 6, only the first and second har- 
monic lines arc detected. During our fitting, we found that 
the uncertainties of third cyclotron lines are very large. So 
we fixed the width at what we got when this parameter is 
freely fitted, which is 1.0, in order to increase the confidence 
of the fitting. The position of the third line is always around 
35 keV, which is not influenced by the fixed width. 

In addition, in some spectra obtained by IBIS and JEM- 
X, the iron line at 6.4 keV is too weak to convince us to add 
this feature into the spectra. Previous authors have found 
that the iron emission line at 6.4 keV is not ncglectablc in the 
spectral fittings (Tammura et al. 1992, Nakajima et al. 2006; 
Mihara et al. 2004 and Tsygankov et al. 2007). Here we have 
adopt the iron lino at 6.4 koV with fixed width at 0.2 kcV 
just as Tsygankov et al. (2007) have done. It doesn't matter 
whether we initially add this additive component into our 
model, because if the iron line is not obvious, then it would 
be too weak to significantly affect the fitting result. 

In Figure 6, the hard X-ray spectra of 4U 0115+63 from 
4-50 keV with RXTE/PCA and RXTE/HEXTE arc pre- 
sented. Here, we select 6 spectra from the well fitted spectra 
of all observational pointings. Indeed all of the spectra ob- 
tained by RXTE have shown prominent iron emission line at 
6.4 keV. Since data of RXTE only covers the energy range 
below 50 keV and spectral data points are not so signif- 
icantly above 30 keV, we only detected the cyclotron ab- 
sorption line features at fundamental line energy around 10 
-15 keV and the second harmonic energy at 22 keV. 



4 DISCUSSION 

4.1 Spin- up Torque Versus Luminosity 

This source is observed continuously during the giant out- 
burst in 2008, which enable us to study the spin properties 
and accretion physics in details. In §. 3.1, we have fitted the 
data points to derive the spin period, and its derivative and 
orbital parameters. We found that the neutron star still un- 
dergone the spin- up process during the 2008 outburst. In the 
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Table 5. Spectral Parameters of 4U 0115+63 with INTEGRAL 



Model Parameters Rev 664 Rev 667 Rev 668 Rev 669 Rev 670 Rev 673 





_ 


1.75 ±1.23 


13.49 ±2.67 


_ 


_ 


_ 


Photo Index 


0.65 ±0.02 


0.15 ± 0.32 


0.20 ± 0.36 


0.65 ±0.02 


0.62 ±0.02 


0.69 ±0.03 


Ecut (keV) 


9.17±0.14 


9.13 ±0.30 


8.23 ±0.77 


9.16 ±0.14 


8.41 ±0.17 


8.42 ±0.23 


Efold (keV) 


12.48 ±0.75 


11.71 ±1,24 


13.13 ±2.51 


12.53 ±0.72 


11.43 ±0.39 


10.37 ± 0.96 


Ecycl,l (keV) 


16.86 ±0.63 


9.53 ±0.62 


10.03 ±1.27 


16.54 ±0.60 


14.65 ±0.47 


15.62 ±0.73 


F-test Probability 


3.87x10-*-* 


1.32xl0-*"'9 


5.343x10-5"' 


6.249x10"™ 


7.566x10-*^ 


1.336 xlO-*'^ 


Depthcyci,! 


0.57 ±0.14 


0.66 ±0.15 


0.73 ±0.31 


0.52 ±0.12 


0.33 ±0.11 


0.58 ±0.18 


Widthcyci,! (keV) 


5.76 ±1.18 


5.25 ±0.68 


5.23±1.45 


5.41 ±1.15 


4.44 ±1.20 


5.11 ±1.52 


Ecycl,2 (keV) 


22.93 ±0.54 


20.57 ±0.31 


20.46 ±0.46 


22.62 ±0.59 


20.97 ±0.99 


22.62 ±0.90 


F-test Probability 


0.0019 


1.09x10-24 


2.377x10-3^ 


0.0017 


0.0037 


5.914 xlO-® 


Depthcyci,2 


0.62 ±0.17 


1.07 ±0.12 


1.09 ±0.20 


1.12 ±0.04 


0.74 ±0.01 


1.10±0.21 


Widthcyci,2 (keV) 


3.42±1.02 


7.51 ±0.59 


6.90 ±0.82 


3.87 ±0.88 


6.14 ±0.85 


4.10±1.25 


Ecycl,3 (keV) 


34.26 ±0.49 


36.69 ±0.33 


35.30 ±0.11 


34.32 ±0.48 


35.09 ±0.76 


35.40 ±1.17 


F-test Probability 


0.0057 


5.213x10-* 


3.421x10"^ 


0.0004 


0.0008 


0.0425 


Depthcyci,3 


0.57 ±0.14 


0.58 ±0.09 


0.38 ±0.11 


0.56 ±0.14 


0.50 ±0.16 


0.98 ± 0.38 


Widthcycl^3 (keV) 


1.00 (fixed) 


1.00 (fixed) 


l.OO(fixed) 


1.00 (fixed) 


1.00 (fixed) 


l.OO(fixed) 


Ecycl,4 (keV) 


43.97 ±2.07 


46.50 ±2.16 


43.74 ±5.58 


43.96 ±1.64 


43.56 ±2.07 


46.05 ±1.58 


F-test Probability 


0.0001 


7.521x10-2" 


9.667x10-22 


0.0001 


0.0005 


0.0297 


Depthcyci,4 


0.73 ±0.13 


1.11 ±0.23 


1.38 ±0.43 


0.74 ±0.14 


0.41 ±0.01 


1.35 ±0.35 


Widthi.yai,4, {keV) 


13.90 ±4.26 


20.48 ±4.38 


24.87 ±7.50 


14.03 ±4.33 


8.49 ±4.54 


9.96 ±5.79 


Ep^ (keV)' (fixed) 


6.4 


6.4 


6.4 


6.4 


6.4 


6.4 


Widthpe (keV) (fixed) 


0.2 


0.2 


0.2 


0.2 


0.2 


0.2 


flux(10~®er(;s cm~^ s~^) 


(1) 5.96 


12.89 


11.05 


9.97 


8.24 


3.04 


(d.o.f) 


2.16(178) 


2.19(182) 


1.63(183) 


2.16(178) 


2.23(181) 


1.69(157) 



<1) range: 5 ~ 50 kcV 



Table 6. Spectral Parameters of 4U 0115+63 with RXTE 



Model Parameters 93032010101 93032010206 93032010301 93032010406 93032010500 93032060201 



Photo Index 


0.32 ± 0.01 


0.40 ± 0.01 


0.16 ± 0.09 


0.29 ± 0.01 


0.31 ± 0.01 


0.56 ± 0.04 


Ecut (keV) 


8.97 ± 0.05 


8.42 ± 0.04 


8.39 ± 0.08 


8.20 ± 0.06 


8.30 ± 0.07 


9.09 ± 0.08 


EfoldikeV) 


8.01 ± 0.18 


8.47 ± 0.07 


8.06 ± 0.23 


8.42 ± 0.14 


8.29 ± 0.19 


8.93 ± 0.54 


Ecycl,i (keV) 


15.42 ± 0.23 


15.23 ± 0.11 


10.04 ± 0.79 


14.34 ± 0.20 


14.66 ± 0.25 


15.55 ± 0.20 


F-test Probability 


1.569x10-20 


6.188x10-1"' 


2.560x10-2" 


5.313x10-3" 


1.31x10-33 


1.926x10-1'' 


Depthcyci,! 


0.30 ± 0.06 


0.43 ± 0.02 


0.41 ± 0.04 


0.32 ± 0.05 


0.28 ± 0.09 


0.16 ± 0.09 


Widthcyci,! (keV) 


2.84 ± 0.46 


4.97 ± 0.32 


7.25 ± 0.76 


3.98 ± 0.45 


3.67 ± 0.72 


2.22 ± 1.13 


Ecycl,2 (keV) 


20.27 ± 0.63 


22.31 ± 0.39 


21.70 ± 0.21 


20.20 ± 0.70 


19.82 ± 1.20 


19.98 ± 1.55 


F-test Probability 


0.0001 


1.333x10-5 


1.966x10-5 


4.049x10-" 


0.002 


2.261x10-" 


Depthcyci,2 


0.64 ± 0.05 


0.24 ± 0.03 


0.39 ± 0.02 


0.53 ± 0.04 


0.52 ± 0.07 


0.79 ± 0.09 


Widthcyci,2 (keV) 


5.29 ± 0.73 


3.96 ± 0.69 


4.59 ± 0.47 


6.64 ± 0.75 


6.62 ± 1.08 


8.43 ± 1.87 


Epe (keV) (fixed) 


6.4 


6.4 


6.4 


6.4 


6.4 


6.4 


Widthpe {keV) (fixed) 


0.2 


0.2 


0.2 


0.2 


0.2 


0.2 


fiux(10— ^ergs cm— 2^- ^) 


(1) 7.06 


15.23 


18.60 


10.46 


9.68 


2.19 


(d.o.f) 


1.65 (57) 


1.98 (57) 


1.99 (57) 


1.58 (57) 


1.41 (71) 


1.06 (57) 



(1) range: 5 ~ 50 kcV 



fittings, we have assume a constant spin-up torque P. How- 
ever, in tiie standard accretion model, tlie spin-up torque 
should correlate to its accretion luminosity. If we could build 
a simplest picture of the disk accretion, material is captured 
to magnetic field lines at and transported to the mag- 
netic poles. We assume all the angular momentum carried 
by the material is transferred to the neutron star, for a Ke- 
plerion disk, the neutron star would have a spin up toque 



(5) 



A neutron star with moment of inertia would subject 
to the torque in Equation ((5]) , when enduring outburst with 
high luminosity, its spin up rate is given by 



= (GM, 



^3/7^2/7^^6/7^ 



(6) 



Assuming the exhibited energy is from accretion, then 
all the gravitational potential energy of the accreted ma- 
terial would transport to the neutron star's surface, with 
luminosity 



GMM^/r 



(7) 



Here is the radius of neutron star. Hence we could 
deduce 



r I 

aL 



6/7 



27r/4GM,)3/7 



(8) 



If we adopt the parameter of neutron star as follows, 
radius ~ 10km, mass Mj, = IAMq, moment iner- 
tia ~ 2Ma;r2/5. The magnetic field of the neutron 
star could be estimated through cyclotron absorption en- 
ergy we've observed. The difference between the first and 
second harmonic cyclotron energy is 10 keV, if we adopt 
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Figure 5. The hard X-ray spectra of 4U 0115+63 obtained by INTEGRAL/ JEM-X (the red line) and INTEGRAL/IBIS (the black 
line) during the observation in revolution 664, 667, 668, 669, 670, 673, the parameters are listed in Table 5 



AE ^ 10fcey"_Bi2(1.2/l +z),z = G^M^/r^c^, hence z~0.2 
and B~ 10^^ G. Then we could get the value of q ~ 
3.348 X 10^"''. 

Additionally, we could estimate the value of a using 
the observed data during the outburst in 2008. We can inte- 
grate Equation (|8]) into Equation ([l]), here the P = = 
-P^u = aP'^L^'\ then we fit the function with the data 
when the luminosity is declining from the peak. The mea- 
sured a = (1.918 ± 0.339) x 10""" is still similar to though 
is a little lower than the theoretical value, suggesting that 
the accretion model described above could well predict the 
spin-up processes of this neutron star system. It should also 



be pointed out that there still exist some uncertainties in es- 
timating a. In observations, the distance of the source and 
the surface magnetic field of the neutron star still have un- 
certainties. In theory, the Eq. (7) assumes that the whole 
accretion power has been transported to X-ray luminosity, 
but in realty only part of the accretion power could trans- 
port to X-ray luminosity. Anyway, using the observed data 
during the giant burst of 4U 0115+63, we have firstly con- 
firmed the relation of spin-up torque versus luminosity in 
this source. 
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Figure 6. The hard X-ray spectra of 4U 0115+63 in the range of 4 - 50 keV obtained by PCA (the black line) and HEXTE (the red 
line) during the observation in pointings 93032010101, 93032010206. 93032010402, 93032010406, 93032010500, and 93032160201. 



4.2 Apsidal Motion 

In this paper, we derived tiie apsidal motion of the binary 
system by comparing different measurements. The apsidal 
motion can probe the gravitational field and the internal 
structure of the companion star. The methods and equa- 
tions (9 ~ 11) follow other authors, e.g. Kelley et al. (1981) 
who discussed the apsidal motion in details, and Raichur & 
Paul (2010) measured three sources' apsidal motion. They 
assume the spin angular momentum of the companion star 
is parallel to the orbital angular momentum, uj is related to 
the apsidal motion constant and the rotational velocity of 
the companion star as (Kelley et al. 1981) 



^ = |?(-f[15g/(e)+ 0^(1 + 9)5(6) 
with 



(9) 



/(e) = (l + fe' + ^e')(l-e^)-', 
<7(e) = (1 - e^)-^ 



(10) 
(11) 



Here g is the ratio of the mass of neutron star to the 
mass of its companion, Rc is the radius of the companion 
star, fl is the ratio of rotational angular velocity of the com- 
panion to the orbital frequency {2n /Porb), o. is the full semi- 
major axis of the orbit. The apsidal motion constant k de- 
scribes the companion's internal matter density distribution. 
We adopt the parameter suggested by Vacca et al (1996) and 
Negueruela & Okazaki (2001), The companion, mass of V635 
Cas is about 19Mq , with radius at 8Rq , i — 43° , rotational 
velocity i^sin i = 290 ± 50 km s-^, hence we could deduce 
k = 0.0025 ± 0.0003. However, there are 70% uncertainty 
in mass and 30% uncertainty in radius of the companion 



Cyclotron Resonance Energies and Orbital Elements of 0115+63 11 



star. If we include these uncertainties, it would lead to a 
larger error bar of 0.0041. Anyway the apsidal motion of the 
binary can put the constraint on the mass-radius relation 
of the companion star. We derive lj ~ 0.048° ± 0.003° ?/r~^ 
which is consistent with the value by Tamura et al (1992) 
and Raichur & Paul (2010). 

4.3 The Resonance Energy and Other Spectral 
Parameters 

In §3, we have obtained all the fitted spectral parameters 
of the X-ray binary 4U 0115-^63 from both INTEGRAL 
and RXTE's data sets covering the outburst of 4U 0115-1-63 
in 2008 April - May, when this source's X-ray luminosity 
changed from (1.28 — 10.90) x lO^^ergss"^ in the range of 
5-50 keV, assuming the distance is 7 kpc. In this section, 
the possible correlations between cyclotron absorption line 
parameters and continuum spectral parameters will be ap- 
proached. Specially, the variations of the fundamental line 
energy during the giant outburst and the possible physical 
origins will be briefly discussed. 

It is still hard to tell which factors could determine the 
cyclotron feature or the continuous spectra. In order to get 
an overview of all the related parameters in spectral fittings, 
we first present some parameters' variations with observed 
time during this giant outburst. 

In Figure 7 we showed the variations of the fundamen- 
tal line energy Eai and the photo index T with observed 
time. The luminosity changes during the outburst are also 
displayed (bottom panel of Fig. 7), where we can clearly see 
the rising phase, the peak phase and deceasing phase of the 
whole outburst. In the plots, we have shown both the spec- 
tral results by INTEGRAL and RXTE. The mean photon 
index V is different for the two instruments during the ris- 
ing and decreasing phases: F ~ 0.6 for INTEGRAL data and 
r ~ 0.4 for RXTE. This effect may be due to the possible dif- 
ferent observing off-angle of the two satellite when observing 
the outburst and there are more absorption lines detected 
in INTEGRAL'S spectrum than RXTE, which would also 
make INTEGRAL'S luminosity a little lower than RXTE's. 

Anyway, the interesting correlation and transition are 
discovered in our observed results. F and Eai vary nearly si- 
multaneously, larger values (softer spectrum) of F at Eai ~ 
15 keV while lower values (harder spectrum) of F at Eai ~ 
10 keV. This correlation is valid for both INTEGRAL and 
RXTE data sets (also see Fig. 8). More importantly, the sud- 
den changes of both parameters occurred during the peak 
phase of the outburst, indicating that there would be a tran- 
sition of this process in the rising phase to the peak phase 
and in the peak to decreasing phase. In the rising phase and 
decreasing phase, Eai ~ 15 keV with a softer spectrum; and 
in the peak phase, Eai ~ 10 keV with a harder spectrum. 
We also noticed that in the spectral fitting of INTEGRAL 
data from 3 - 100 keV, only during the peak phase, the 
significant column density absorption is detected (see Table 
5). 

This sudden changes of spectral property and funda- 
mental line energy suggest that some physics results in this 
transition. The possible correlation between X-ray luminos- 
ity and Eai have been pointed out before. We also dis- 
play the possible relation between X-ray luminosity and 
Eai in Fig. 9. As we mentioned previously, the derived X- 
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Figure 7. The evolution of the fundamental resonance energy 
and photo index with time by both INTEGRAL (blank squares) 
and RXTE (solid squares) observations. The sudden changes of 
both parameters are observed from the transition in the rising 
phase to the peak phase and in the peak to decreasing phase. 



ray luminosity from INTEGRAL data is generally lower 
than that from RXTE data (by ~ 20%). For each instru- 
ment's results, we can see the possible correlation: higher 
luminosity with lower values of Eai- But the sudden tran- 
sition of the Eai changes are clearly found rather than 
a real correlation function. For INTEGRAL'S data sets, 
when Lx > 6.76 x lO'^^ergs s~^, we detect the familiar fun- 
damental CRSF at Eai ~ 10 keV, below this luminosity 
level, the fundamental line energy changes to be Eai ~ 15 
keV; while as for RXTE's data sets, the criterion is around 
Lx > 8.5 X lO^^ergss"^. For RXTE data points, we also 
noticed that for the similar X-ray luminosity ranges (from 
(8.4 — 9) X 10^^ ergs s~^), both two fundamental line ener- 
gies at ~ 10 keV and ~ 15 keV are detected, suggesting that 
the X-ray luminosity may be not the unique factor in deter- 
mining the fundamental line energy variations. This transi- 
tion is not complete yet because we do not get enough data 
points when the fundamental line energy distributes from 
Eai ~ 12 — 14 keV during the sudden transition phases. 

Burnard et al. (1991) studied this luminosity dependent 
relationship, and assumed that the CRSF is formed above 
the surface of neutron star at height hr in the accretion 
column. Based on the dipole law of the magnetic field, the 
fundamental cyclotron absorption should follow the equa- 
tion given below: 



12 Jun Li, Wei Wang and Yongheng Zhao 




1 ' i I I 1 * 1 I I 



^ 



I I I 



Lx(10 erg/s) 



Figure 8. The relationship between photo index and fundamen- 
tal cyclotron energies (data points obtained from both INTE- 
GRAL and RXTE). 



Figure 10. The height, converted from equation (13) versus lu- 
minosity. The dashed line represents the column top height, cal- 
culated by equation (14). 




Lx (10 erg/s) 



Figure 9. The relationship between X-ray luminosity and fun- 
damental cyclotron energies (data points obtained from both IN- 
TEGRAL and RXTE). 



Eai oc {Rns + K) -^l + Zg)- 



(12) 



where Rns stands for the radius of the neutron star, Zg 
denotes for the gravitational redshift. Here we could treat 
the redshift item as a coefficient, for simplicity, we set Zg — 0. 
Then we could describe the influence of height on resonance 
energy in a convenient way. 

hr 



Rns 



i—y 



1/3 



Burnard et al. (1991) have given their estimation: 



reff TT 



-Rns- 



(13) 



(14) 



^"Edd ^® Eddington luminosity along the magnetic field, 
identical to the conventional Eddington luminosity for a neu- 



tron star with a mass of 1.4M©, L 



2.0 X W^^ergss- 



H± is the ratio of the Thomson cross section to the 
Rosseland-averaged electron scattering cross section for ra- 
diation flows across the magnetic fields. Ifere we adopt H± 
= 1.23 deduced by Mihara et al. (2004). 

So we could employ the Eai into the above equations, 
then estimate hr/Rns as a function of X-ray luminosity 



from our observations by INTEGRAL and RXTE. The re- 
sults are shown in Figure 10. Here we assume -Eo = 18 keV, 
which is the fundamental resonance energy at the surface of 
the neutron star, taken to be the possible maximum value 
of our observations. We deduce all the hr/ Rns according to 
equation (13). However we did not find the similar results 
done by Mihara et al. (2004). Even we have treated data 
sets from INTEGRAL and RXTE respectively, our results 
could not confirm this theoretical model. 

Instead of this linear relation ship suggested as equation 
(14), we find a phase transition similar to F and Eai- The 
trigger of the sudden transition may closely related to the 
luminosity of the source. As the luminosity increases above 
a certain level, condition might change suddenly instead of 
what equation (14) has suggested. Besides, as shown in Fig- 
ure 10, the hr/RNS ratio drops suddenly when it decreases 
to under 0.1 level, this trend could also be seen in the work 
done by Nakajima et al. (2006). So the mechanism of the 
fundamental cyclotron absorption may be affected by other 
unknown physical factors, or the magnetic field at the sur- 
face of neutron star might not be consistent with the dipole 
assumption. 

We also study the ratios EaijEai and hr/RNS ver- 
sus time of the outburst, as shown in Figure 11. Since Ea\ 
changes significantly during the outburst, while Ea2 is rela- 
tively stable around 22 keV in the whole outburst, the vari- 
ation of the ratio Eai/Ea2 shows the similar behavior to 
the evolution of Eai- in the rising phase, the ratio is about 
1.3, in peak phase the ratio switches to be around 2, and in 
decreasing phase, the ratio returns to be near 1.3. The tran- 
sition may be related to the X-ray luminosity level, deserved 
more detailed studies. 

Coburn et al. (2002) found that there exists a positive 
correlation between absorption depth and W/Ea in phase 
averaged spectra of the X-ray pulsar. This indicates that this 
source has a tall cylindrical shape accretion column rather 
than a fiat coin shape, suggested by Kreykenbohm et al. 
(2004). In Figure 12 we presented the ratio of the width 
to energy of both Eai and Ea2 versus the depth. The posi- 
tive correlation between W/Ea and depth is found for both 
two lines, though large scattering and uncertainties still ex- 
ist. This result, along with what Kreykenbohm et al. (2004) 
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Figure 11. The evolution of Ea2/ Ea\ and hr/ Rns with time of 
the outburst. The circles represent the luminosity of the outburst. 
We could see the sudden changes at the transition from rising 
phase to peak phase in MJD 54554, and transition from peak 
phase to decreasing phase in MJD 54562. 
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Figure 12. W/Ea of fundamental and second harmonic reso- 
nance energy of 4U0115-I-63, varies with their energy values, re- 
spectively. 



and Nakajima et al. (2006) have obtained, strengthens the 
cyUndrical column geometry assumption. 

We have fitted aU the continuum spectra with a power- 
law model plus a high energy cutoff. The photon index vari- 
ations have been plotted before. Finally, we presented the 
variations of other two spectral parameters Ecut and Efoid 
with time during the outbursts in Fig. 13. Most values of 
Ecut do not change too much during the whole outburst, 
it stays around 8.5 keV for both RXTE and INTEGRAL 
observations, which is consistent with the results by Tsy- 
gankov et al (2007) and Makishima et al (1999). For Efoid, 
INTEGRAL results show a higher values (~ 12 keV) than 
that (~ 8 keV) obtained by RXTE, which would be due to 
the wider energy range by INTEGRAL which covers higher 
ranges up to 100 keV. 



a> a 



LU 8 
7 

14 



0) 12 



10 



12 

^ 10 

D) 8 
0) 

"o ' 

^ 4 

X 

_] 

2 





■ INTEGRAL 
□ RXTE 



■ INTEGRAL 
□ RXTE 



n 5 



O 



54540 54550 54560 54570 54580 

MJD (day) 



Figure 13. The evolution of cut-off energy and fold energy with 
time during the outburst (data points from both INTEGRAL and 
RXTE). 



5 SUMMARY AND CONCLUSION 

The giant outburst in neutron star high mass X-ray binaries 
provides us a good chance to study the accretion physics and 
interaction between magnetosphere and accretion materials. 
We observed the giant outburst of 4U 0115-1-63 during April 
- May 2008 using both INTEGRAL and RXTE data, per- 
forming both timing and spectral analysis in details. The 
spin period is determined to be 3.61430 ± 0.00003 s at MJD 
54566, and we confirm that this source is enduring spin up 
state at rate of Apin = (-7.24 ± 0.03) x lO'*' s d~^ dur- 
ing the outburst. In addition, we determined the angle of 
periapsis of the orbit oj = 48.67° ± 0.04° during the 2008 
outburst, as we compare with previous works done by other 
authors, we get lu = 0.048° ± 0.003° yr"\ We also noticed 
the pulse profile's merging behavior during the outburst. 

Furthermore, spectral properties of 4U 0115-1-63 dur- 
ing the different phases of the outburst are studied. All 
spectra from both INTEGRAL and RXTE are fitted with 
the power-law model plus a high energy cutoff and added 
by two components, cyclotron resonant absorption features 
and a gaussian line at 6.4 keV. We have searched for the 
cyclotron resonant absorption features up to 5th harmonic 
cyclotron line with INTEGRAL data. The fundamental line 
energy varies during the outburst at 10 keV or 15 keV re- 
spectively; the second harmonic energy at 22 keV; the third 
at ~ 34 keV; 4th at ~ 45 keV. The fifth harmonic energy 
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is not significantly detected. ISence we have confirmed that 
4U 0115+63 shows the four significant cyclotron absorption 
lines features. 

The present observations with INTEGRAL and RTXE 
provide us an opportunity to study variations of all spec- 
tral parameters in this source with time and resonance en- 
ergy during the 2008 giant outburst. A positive relation- 
ship between photo index and fundamental cyclotron energy 
is found (soo Figure 8). The evolutions of fundamental cy- 
clotron energy and spectral properties with time during the 
outburst show the sudden transitions: at the rising phase, 
the source shows the soft spectrum with Eai ~ 15 keV; 
in the peak phase, Eai ~ 10 kcV with a harder spectrum, 
and high column density; then into the decreasing phase, 
returning to the soft spectrum with Eai ~ 15 keV. The 
luminosity dependence of resonance fundamental energy is 
also indicated in our observations, but in the similar lumi- 
nosity range of ~ (8.4 — 9) x 10^^ ergs s~^, both two fun- 
damental energies at ~ 10 and 15 keV are detected. So the 
physical parameters other than the single X-ray luminos- 
ity may play a role in determining the spectral properties 
during the outburst. Instead of a positive relationship of 
hr/RNS and luminosity as expected, we also found a transi- 
tion which makes a sudden change for the ratio (below ~ 0.1 
or around 0.2), similar to fundamental line energy, which in- 
dicates that the environment near the surface of neutron star 
is much more complicated than we originally thought. This 
transition process during the outburst may require further 
studies both in observational and theoretical work. The X- 
ray luminosity should be a important parameter, when it 
reaches a certain critical value, sudden changes of accre- 
tion and radiation could occur or other physical conditions 
would change. Another positive correlation between depth 
and W/Ea ratio is consistent with previous work, which sup- 
ported the tall cylindrical geometry assumption of the accre- 
tion column. The distinct spectral properties (different pho- 
ton indices and fundamental line energies) between the peak 
phase and the rising/decreasing phases suggested different 
accretion physics and environments in the peak luminosity 
ranges. 
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